Modeling inelastic phonon scattering in atomic- and molecular-wire junctions 
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Computationally inexpensive approximations describing electron-phonon scattering in molecular- 
scale conductors are derived from the non-equilibrium Green's function method. The accuracy is 
demonstrated with a first principles calculation on an atomic gold wire. Quantitative agreement 
between the full non-equilibrium Green's function calculation and the newly derived expressions is 
obtained while simplifying the computational burden by several orders of magnitude. In addition, 
analytical models provide intuitive understanding of the conductance including non-equilibrium 
heating and provide a convenient way of parameterizing the physics. This is exemplified by fitting 
the expressions to the experimentally observed conductances through both an atomic gold wire and 
a hydrogen molecule. 
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The rapid evolution in electronics towards smaller and 
faster devices will eventually reach the fundamental level 
set by the atomistic structure of matter. Atomic-size con- 
ductors take this development to the extreme of minia- 
turization Q, and understanding their properties is an 
important problem in the emerging fields of nanoelec- 
tronics and molecular electronics. One relevant aspect 
is the study of the effects caused by atomic vibrations, 
since inelastic scattering of traversing electrons and en- 
ergy dissipation play essential roles for device character- 
istics, working conditions, and stability. Vibrational sig- 
nals can also be used to extract information about the 
detailed microscopic configuration, which usually cannot 
be imaged simultaneously with a transport measurement. 
Inelastic effects have in the recent years been studied in 
a variety of nanoscale systems, e.g., single molecules on 
surfacesprobed with the scanning tunneling microscope 
(STM) 0, molecules in break junctions Q, and metallic 
atomic wires 

Theoretical descriptions of inelastic transport through 
small devices connected to metallic contacts include 
many-body theory in the Coulomb blockade regime 
single-particle first-order perturbation approaches 
i.e., "Fermi's golden rule" (FGR), as well as calcula- 
tions to infinite order based on the self-consistent Born 
approximation (SCBA) combined with non-equilibrium 
Green's functions (NEGF) [111 OH 0ur work is based 
on the SCBA, which in contrast to FGR takes the many- 
particle nature of the problem into account. However, 
the SCBA method is computationally very demanding 
especially when used in combination with first principles 
electronic structure methods. Moreover, the SCBA does 
not yield simple formulas which can be used to extract 
information from experimental data. 

In this paper we develop methods which vastly simplify 
the SCBA approach. The main results are analytical for- 



mulas for the current and power derived from a lowest 
order expansion (LOE) of the SCBA expressions. In par- 
ticular, we show how first principles SCBA calculations 
on atomic gold wires can be accurately described by the 
LOE with minimal computational effort. Moreover, we 
derive compact analytical expressions using two simple 
models. These latter models are able to fit both the the- 
oretical SCBA results as well as experiments using the 
electron-hole damping rate of the phonon as the central 
parameter 

Phonon scattering is included in the SCBA method as 
self-energies to the electronic description. We use the un- 
damped phonon Green's functions to express these self- 
energies in the device subspace as 
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Here, M A is the electron-phonon coupling matrix for 
phonon mode A occupied by n\ phonons with energy 
hu>\. The lesser /greater self-energy matrices £> h are 
given by two terms corresponding to absorption/emission 
of phonon quanta. We furthermore assume that 
these self-energies can be used in non-equilibrium with 
a bias dependent phonon occupation number n\(V). 
The retarded self-energy can then be obtained from 
the greater/lesser parts using the Hilbert transform 
(H {/(£')} (E) = l/nVj f(E')/(E E>) dE'). 
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1 The polaron term |l7l in the retarded self-energy in Eq. J2J has 

been neglected since it gives no "signal" at the phonon energy. 
However, it gives rise to two additional terms in the expression 
for the current (Eq. proportional to V and V 2 and does not 
contribute to the power. 
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The computational difficulty of solving the SCBA 
equations stems from the coupling of Green's functions in 
energy. Calculations usually involve a numerical energy 
grid which has to be fine enough to resolve the low tem- 
perature structure of the Fermi function, while at the 
same time span a large energy range to cover phonon- 
energies, applied bias, and allow an accurate computa- 
tion of the Hilbert transform which is nonlocal in energy. 
The current and power are then computed as integrals 
over this energy grid 0, 0, 0] . 

These difficulties can be overcome if (i) the electron- 
phonon coupling is weak, i.e., the probability for multi- 
phonon processes is low, and (ii) the density of states 
(DOS) of the contacts and the device is slowly varying 



over a few phonon-energies around the Fermi energy Ep, 
i.e., in the notation used below, G r (E) rs G r (Ep ) and 
Ti y2 (E) w Ti^(Ep). These approximations are valid 
for systems where (i) the electron spends a short time 
compared to the phonon scattering time in the device and 
(ii) the closest resonance energy (E Tes ) is either far away 
from the Fermi energy (\E Tes — Ep\ ^> T, eV and Huj) or 
the broadening by the contacts is large (r » eV, Huj and 
[E m — Ep\). The expressions for the current and power 
can then be expanded to lowest order (second) 
in the electron-phonon coupling and the integration over 
energy performed analytically. The power dissipated into 
the phonon system p LOE can, after lengthy derivations, 
be written: 



pL ° E = ^2^-^-Mfiux)-nx) Tr[M A AM A A]+P(y,^ A ,T)Tr[M A Gr 1 G t M A Gr 2 Gt] , (3) 

A 

= Huj (cosh (H) - 1) coth ( Huj - eV sinh {%) 
nh cosh (^)_ cosh (|V) 



where ub is the Bose-Einstein distribution, which ap- 
pears naturally from the integration of the Fermi func- 
tions of the electrons in the contacts. Here, G = G r (Ep) : 
Ti i2 = Ti, 2 {Ef), and A = i(G — G^) are the non- 
interacting, i.e., without phonon interactions, retarded 
Green's function, the broadening by the contacts, and 
spectral function at Ep, respectively. 

From Eq. © we see that the power can be decom- 
posed into terms corresponding to the individual phonon 
modes. We also note that the first term describes the 



power balance between the electron and phonon sys- 
tems (at zero bias) with an electron-hole damping rate 
7g h = cj A /7rTr [M A AM A Aj and is in fact equivalent to 
the FGR expression [ill Il4| . The second term is even in 
bias and gives the phonon absorption/emission at non- 
equilibrium; it is negligible at low bias (eV <C Huj), turns 
on at the phonon energy and becomes linear in voltage 
at high bias (eV 3> Huj). 

Using the same approximations, the current through 
the device J LOE is given by [l5| : 
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where np is the Fermi function, the bias is defined via, eV = /i 2 — fii, and the conductance quantum Go = e 2 /irH 
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FIG. 1: Universal functions (Eqs. JSJ and iJTJ) giving the 
phonon contribution to the current. The differential conduc- 
tance dl/dV and the second derivative signals is shown for 
one phonon mode with the bias in units of the phonon energy 
at a temperature kT — 0.025 hu). For the symmetric term, 
the FWHM of the second derivative peak is approximately 
5.4fcT:i8|. 
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FIG. 2: Comparison between the SCBA results and the LOE 
expressions (Eq. J^J) without heating a) and with heating b) 
(7d = 0) at T = 4.2 K. The parameters for the ABL model 
(Eq. il'21 ) were extracted directly from the DFT calculations, 
7eh = 5.4 x 10 10 s _1 and hu; = 13.4 meV. 



appears naturally. In contrast to the 1st order Born ap- 
proximation, these expressions are current conserving like 
SCBA. 

The current expression retains the structure of the 
Landauer expression (first term of Eq. (J5J) and gives cor- 
rection terms for each phonon mode. The phonon terms 
can in turn be divided into a "symmetric" term Z s y m 
where the differential conductance dl/dV is even in bias, 
and an "asymmetric" term containing the Hilbcrt trans- 
form Z Asym yielding an odd contribution. We note the 
simple factorization into terms depending on the elec- 
tronic structure at Ep and universal functions 2 Sym and 
-jAsym wn j c h yield the line-shape of the inelastic signals in 
the I—V, see Fig. ^ Whether the conductance increases 
or decreases due to phonon scattering depends on the 
sign of the traces in Eq. JSJ and will be discussed further 
below. Examination of the "asymmetric" term in Eq. (JSJ 
shows that it is zero for symmetric systems. Although ex- 
perimentally measured conductances contain asymmetric 
signals, the size of these signal is usually small in the pub- 
lished curves. At present it is unclear if they are caused 
by phonons or other effects. 

As we have shown previously heating of the phonon 
system should be considered [9j which makes the num- 
ber of phonons n\ bias dependent. The simplest way to 
include non-equilibrium heating is to write down a rate 
equation, including an external damping rate 7^ of the 
phonons: 

pLOE 

hx = -± + 7d (n B (thu x ) - n x ) , (8) 



phonon modes 2 . The steady state occupation nx is eas- 
ily found. Substituting the result into Eqs. 10)- gives a 
computationally simple but powerful formula for the cur- 
rent through the device including heating of the phonon 
system. 

To judge the accuracy of the LOE approach, we com- 
pare the LOE results to the full SCBA solution for a 
four atom gold wire, see Fig. [21 The SCBA calcula- 
tion was performed as described previously |9j , where the 
Hamiltonian, phonon modes, and electron-phonon inter- 
action were obtained from density functional calculations 
(DFT). The excellent agreement between the full SCBA 
and the LOE expression can be understood by noting 
that the DOS of a gold wire is slowly changing over an 
energy range much greater than the phonon energies. In 
addition, the electrons only spend a small time in the wire 
compared to the electron-hole damping rate. Impor- 
tantly, the LOE conductance calculations were performed 
in less than a minute on a regular PC, compared to sev- 
eral hours for the SCBA calculations. The LOE approach 
thus opens up the possibility to study inelastic scatter- 
ing with first principles methods for large systems, e.g., 
organic molecules. 

To gain further insight into the expressions presented 
above, we consider a single electronic site with symmetric 
contacts r = Ti = T2 coupled to one phonon mode. 
Introducing the transmission probability r = |G| 2 r 2 and 
the electron-hole damping rate 7 c h = 4(cj/7t) M 2 t 2 /T 2 , 



where pL° E is the power dissipated into the individual 



2 For weak electron-phonon interaction, the division of power into 
the individual phonon modes is straightforward from Eq. i fl 
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FIG. 3: a) Single level model (Eqs. <l9TTuT l fitted to the 
experimentally measured conductance through a Deuterium 
molecule |T^I . The parameters used for the fit are hu> = 
50meV, r = 0.9825, 7 ch _ 1.1 x 10 12 s -1 , and T = 17K. b) 
The ABL model (Eqs. I|llll2|l i fitted to the measured con- 
ductance through an atomic gold wire (experimental data 
from Ref. 4\). The fit reveals the following parameters, 
hu = 13.8 meV, T = 10 K, 7 ch = 12xl0 10 s" 1 , and j d = 37 ch . 



we obtain: 
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rLOE 



7 ch frw (n B (hw) - n) + ^ J—^, (9) 

rV + e lch — — X Sym . 10 

4 erko 



We note that, from the term 1 — 2t in Eq. (|10f) . the con- 
ductance will increase due to phonon scattering for low 
conductance systems (r < 1/2) and decrease for highly 
conducting systems (t > 1/2). The LOE approach di- 
rectly provides the sign of the conductance change in 
contrast to FGR approaches where this requires careful 
considerations [flul- 

The conductance through a single hydrogen molecule 
has been measured using a platinum break junction setup 
0, 0| . Because the elastic current is carried through a 
single molecular orbital ,16], the single level model fits 
the experiment very well, see Fig. The best fit is ob- 
tained using a negligible external damping of the phonon 
mode (7d -C 7ch) which can be understood physically 
from the mass difference between the hydrogen molecule 
and the platinum atoms of the break junction. Wc also 
note that both the size of the conductance step and the 
conductance slope (caused by heating) is fitted with only 
one parameter, the electron-hole damping rate 7 e h- 

The electronic structure of atomic gold chains are qual- 
itatively different from the one level model. However, it 
is relatively straightforward to derive an alternating bond 
length (ABL) model. Inserting the electron-phonon ma- 
trix for an ABL phonon mode and using the Green's 
function for a half filled perfectly transmitting 1-D chain 



^ABL = 7ehfiw[n B (fto;) -n] 



7ch 
2 



ftuJ 



V, 



rLOE 
^ABL 



e 2 y e 7ch lSvn 
nh 2 ehu 



(11) 
(12) 



where the only difference to the one-level model is that 
t = 1 (perfect transmission) and a factor of two caused 
by the absence of forward scattering from an ABL mode 
(the one-level model has an equal amount of forward and 
back scattering). The ABL model is shown in Fig. [3 with 
the 7 c h damping rate calculated directly from the DFT 
model. The main difference compared to the SCBA/LOE 
results is the assumption of perfect transmission through 
the chain. Fitting the ABL model to experimental data 
gives the very satisfactory fit shown in Fig. 03. We 
briefly note that the external damping 7d = 3 7 e h is not 
negligible in contrast to the hydrogen case. In this pa- 
per we have used sharp phonon energies, c.f., Eq. J2J. 
However, if the phonon spectral function is known, it is 
possible to introduce broadening directly into Eqs. ®- 
(|7|) from a finite phonon lifetime. 

We have derived simple and accurate approximations 
to describe the effect of phonon scattering on the con- 
ductance through nanoscale conductors. The approxi- 
mate expressions greatly reduce the computational effort, 
compared to solving the SCBA equations. In addition, 
simple models were derived which provide insight and are 
suitable to fit experimental data. 
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